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Abstract 



We explore the tunneling behavior of a quantum particle on a finite 
graph, in the presence of an asymptotically large potential. Surprisingly 
the behavior is governed by the local symmetry of the graph around the 
wells. 

1 Introduction 

Quantum tunneling is the physical phenomenon that a quantum particle 
can get to the other side of an energy barrier. This is of course not possible 
in classical mechanics. Even more counterintuitive, Oskar Klein in [5] 
observed that the higher the barrier is, the higher the probability is that 
the particle crosses the barrier. As the height goes to infinity, the barrier 
becomes invisible to the particle. In this paper we shall investigate to 
what extent and under what conditions does quantum tunneling happen 
in the discrete setting. In particular we shall see whether Klein's paradox 
can be recovered. 

The energy barrier is easiest realized by a potential that has two or 
more local minima, referred to as energy wells. In physics, the motion of a 
quantum particle on a manifold is governed by the Schrodinger equation: 



Here *t is the quantum state of the particle and H = A — V where A is the 
Laplace operator and V is a given potential. Since the Laplace operator 
has a natural discrete analogue on graphs, this makes it possible to define 
and analyize the discrete Schrodinger equation. 

Let G(X, E) be a finite connected graph. For any vertex a; £ A we 
denote its degree by d x . We think of the vertices as the possible positions 
of a quantum particle. At any moment t the state of the particle is given by 
a unit length vector <p t £ C x ,Y^ x ex \ifit(x)\ 2 = 1. The quantity |ip t (a;)| 2 
is interpreted as the probability of the particle being in position x. Let A 
denote the symmetrized Laplace matrix of the graph, that is 
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Given a potential V : X — > [0, oo), the motion of the particle is gov- 
erned by the discrete Schrodinger equation 

- i^ft = Hipt (2) 

where H — A — V is the Hamiltonian. (Here and from now on we will use 
the same notation for a function on the set of vertices and the correspond- 
ing diagonal matrix.) We want to investigate the asymptotic behavior of 
this quantum evolution, so we choose the potential to be V = Q ■ W where 
W is a fixed vector and Q is a real number going to infinity. 

Simon (see [5]) carried out asymptotic analysis of double well poten- 
tials in Euclidean space for rather general potentials. We will consider 
potentials with two (double well) and three (triple well) global maxima. 
We will consider simple potentials meaning that the only vertices where 
the potential is non-zero are the wells themselves, and general potentials. 
It turns out that whether tunneling happens at all is primarily governed 
by the degree of similarity of the neighborhoods of the wells, and the speed 
of tunneling depends on the distance of the wells. More precisely: 

Definition 1.1. Let PR(x,k) denote the probability that the simple 
random walk started at x £ X returns to x at time k. Given two vertices 
x, y € X we say that they are m-cospectral if the probability of return is 
the same for x and y up to time m, that is PR(x, k) — PR(y, k) for every 

< k < m. The cospectrality of x and y is the maximal m for which they 
are m-cospectral. This will be denoted by co(x,y). 

Definition 1.2. Given two states x,y € X, let us start the particle in 
the pure state of x. We define the tunneling coefficient to be 

liminf sup \<p t {y)\ 2 = TC(x,y). (3) 

Q^oo te [o,oo) 

We distinguish three different behaviors: Perfect (asymptotic) tunnel- 
ing happens if TC(x,y) = 1, partial tunneling happens if < TC(x,y) < 

1 and no tunneling takes place if TC(x,y) — 0. 

We say that the tunneling from x to y happens within f(Q) time if 

TC(x,y) = ]hami sup |<My)| 2 . (4) 

Q^°° te[o,/(Q)] 

We will use the Bachmann-Landau notation to talk about tunneling 
times. 

Theorem 1.3. In the case of a simple double-well potential where the 
graph-distance of the two wells x, y is denoted by d — d(x, y) the following 
is true between x and y: 

a) if the two wells are d-cospectral (co(x,y) > d) then there is perfect 
asymptotic tunneling 

b) if co(x, y) = d — 1 then there is partial tunneling 

c) ifco(x,y) < d — 1 then there is no tunneling. 
When there is tunneling, the tunneling time is Q(Q d ~ 1 ). 
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The triple-well setup turns out to be much more complicated. To 
avoid excessive case analysis, here we restrict our attention to sufficiently 
cospectral wells. Let x, y, z £ X be three points such that a = d(x,y) < 
b = d(x, z) < c = d(y, z). Assume that the wells are pairwise c-cospectral, 
and let the potential be W(x) = W(y) = W(z) = 1 and everywhere 
else. For any two wells u, v let us compute 

l-PI-i 

= ^2 n i i \j d *i d *j+i ( 5 ) 

where P : u — > v denotes any path u — xo, xi, . . . xi = v in the graph and 
the length of this path is denoted by \P\ — I. 

Theorem 1.4. With the above notation and assumptions 

a) If a < b < c then there is perfect tunneling between x and y, but no 
tunneling to and from z. The tunneling time is Q(Q a ~ 1 ). 

b) If a = b < c and the particle is started from x then it never tunnels 
perfectly to a single other state, but it does tunnel to a mixed state 
of y and z. This happens on the scale of t ~ Q a_1 . 

c) If a = b < c and the particle is started from y, then it tunnels com- 
pletely to z if and only if c xy = c xz . The tunneling time is B(Q a ~ 1 ). 
Note that in particular, unlike in the first case, the particle can tun- 
nel from y to z despite the fact that d(y,x) < d(y,z). Comparing to 
the double-well case, note that the speed of tunneling between y and 
z increased substantially after introducing the third well at x. 

The precise behavior becomes even more subtle If the wells form an 
equilateral triangle. 

Theorem 1.5. Assume a — b — c. 

a) If c xy 7^ c xz then there is only partial tunneling from y to z. 

b) If c xy = c xz and \Jc yz + 8c xz is irrational, then there is perfect 
tunneling from y to z. The tunneling time f has to satisfy Q a_1 = 

o(f(Q))- 

c) If Cxy = c xz and \Jc yz + Sc xz is rational, then even if there is perfect 
tunneling from y to z, the tunneling time f has to satisfy Q a = 
0(f(Q)). (In particular this is the case if c xy — c xz — c yz .) 

Conjecture 1.6. If c xy = c xz 7^ c yz we believe that there is always 
perfect tunneling from y to z. 

In view of the characterisation of the double-well case, there is a very 
surprising instability phenomenon in the triple-well case. 

Theorem 1.7. It is possible to construct a triple-well scenario where 
TC(y,z) = 4/9, but by modifying the graph arbitrarily far from the 
wells only a little bit (e. g. by adding a single edge) one can achieve 
TC(y,z) = 1. Hence the tunneling coefficient is not predictable from any 
fixed neighborhood of the wells. 

Finally we show that tunneling can happen in the presence of more 
general potentials then the singular ones. In particular perfect tunneling 
happens in the symmetric double-well case. 
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Theorem 1.8. Let G(X,E) be a graph, x,y € X and W : X — > R which 
has its global maxima in x and y. Let us further assume that there is an 
involution of G that takes x to y and that preserves W . Then there is 
perfect asymptotic tunneling between x and y and the tunneling time is of 
the usual order: Q^**)- 1 . 



1.1 Intuitive explanation of tunneling 

It is easy to see that if the inital state of the particle is tp = ipo then the 
solution of the Schrodinger equation is given by 

itH / n \ 

<pt — e (p. (6) 

This evolution (without the potential) is usually referred to as the contin- 
uous time quantum walk. It has been first studied in [2] . For an overview 
the reader is referred to [?]. The main application of quantum walks is 
in quantum computing. Since the quantum walk "moves" faster than the 
classical random walk (see e. g. [T|) it can be used to show the power of 
quantum computing over classical computers. 

To better understand the behavior of the solution we decompose the 
state-space according to eigenstates. Let n = \X\ denote the number of 
vertices of G. The matrix H is symmetric hence it has a real spectrum 
with eigenvalues Ai < A2 < A3 < ■ • • < A n . Let 4>k denote the eigenvector 
corresponding to the eigenvalue Afc. Then e ltH has the same eigenvectors 
but with eigenvalues e ltXk . Thus if we write 



then we get that 



ft = ^c fe e ltAfc V'fc- (8) 



If the potential well is deepest at the vertices x, y £ X with value 
V(x) — V(y) = Q then from the physical analogy one expects the two 
smallest (most negative) eigenvalues to be approximately 1 — Q, and the 
corresponding two normalized eigenfunctions should be close to tpi(x) ~ 

My) « M%) ~ -My) « Vv^ and v>i(z) « M*) « if z e x\{ x , y}. 

Then if the particle starts in x we have 
Hence 

^ t ^^(e lt(Al - A2 Vi+^2) (10) 

so at time t = 7r/(Ai — X2) the state of the particle is approximately given 
by c(ipi — 1P2) so the particle is almost surely in state y. 

In the following sections we shall make precise estimates to be able to 
confirm or reject the predicted behavior of the particle. 
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2 Asymptotics 



In order to make the intuitive explanation work in reality, we have to 
understand asymptotic behavior of various quantities. The initial state of 
our particle will be ip = \ x , the characteristic function of the vertex x G X 
which is one of the vertices where the potential is minimal. Hence first of 
all we will be interested in the eigenvector-decomposition ip = J^™ Ckipk- 
To be able to analyise this, we will first need to understand asymptotics 
of the eigenvalues and eigenspaces of the Hamiltonian. 

We shall always denote the increasing sequence of eigenvalues of H — 
A — Q ■ W by Ai < A2 < ...A n , and the corresponding normalized, 
pairwise orthogonal eigenvectors by ipi, ip2, ■ ■ ■ , ip n - (If there are multiple 
eigenvalues, the choice of ipi might not be unique.) 

2.1 Eigenspaces 

The spectrum of a matrix depends continuosly on the entries. Since 
H/Q = A/Q — W — y —W if Q — y 00 we can observe the following: 

Claim 2.1. If the values of the unsealed potential W are wi < W2 < • • • < 
w n then for all 1 < i < n we have 

lim ^ + Wl = 0. (11) 

Lemma 2.2. For any i the eigenfunction ipi is asymptotically concen- 
trated on those vertices of the graph on which the value of the unsealed 
potential is equal to u>i. That is ipi(x) — > as Q — > 00 unless W(x) — Wi. 

Proof. Denote by the characteristic vector of x G X and write = 
Y^Cktpk- Then ipi(x) = (ipi,Xx) = d- We further have Hxx = I] A fe c fc t/> fc 
and hence 

A X * = (H + QW)x* = H Xx +QW(x) Xx = ]T(A fc +QW(x))c k ip k . (12) 
Dividing both sides by Q and taking the norm-squares we get that 

The left hand side converges to zero, and if W(x) 7^ w» then Ai/Q + 
W(x) converges to a non-zero constant, hence a has to converge to as 
Q goes to infinity, and this finishes the proof. □ 

A level set of the potential W is a maximal subset L C X on which W 
is constant. The preceeding lemma shows that for any such level set there 
will be precisely \L\ eigenfunctions of the Hamiltonian that become more 
and more concentrated on L, and conversely, any function concentrated 
on L will be more and more composed only of these \L\ eigenfunctions. 

The "wells" are the vertices of the level set corresponding to the largest 
value of W (since we are subtracting W this is where the potential is ac- 
tually minimal) . As we are interested in the behavior of a particle started 
from one of the wells, we need to understand the asymptotic behavior 
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of the eigenfunctions corresponding to this level set. We know that they 
become concentrated on the wells and we already understand the first 
order behavior of the corresponding eigenvalues. However to study the 
tunneling phenomenon, we need a subtler description of the eigenfunc- 
tions involved to see whether tunneling actually happens or not. We also 
need a more precise asymptotic on the eigenvalues to be able to estimate 
the tunneling time. 



2.2 Eigenfunctions 

The method we use to understand eigenfunctions comes from the theory of 
random walks. It is well known that for any nonempty subset L C X and 
a fixed function / : L — > R there is a unique extension / : X — > R that is 
harmonic everywhere outside of L. One possible way to do this extension 
is to start a random walk from a vertex x and define f(x) to be the 
expected value of / at the point where the walk first hits L. This method 
can be extended to the case, where we want / to be an eigenfunction of 
the Laplacian with eigenvalue A outside of L. Denote the vertices of the 
random walk by x = vo , vi , . . . , vt where T is the first time when the 
walk enters the subset L. Then define f(x) to be the expected value of 
/(ut)/(1 — A) T . The prescence of the potential makes things even more 
complicated, but the construction below is still inspired by the previous 
argument. 

Let L C X be the level set of the wells. We may assume without loss 
of generality that W\l = 1 and < W^|x\£ < 1. Fix an arbitrary function 
/ : L — y R and extend it to X \ L the following way. Let us consider a 
walk P = (xo,xi, . . . ,xt) of length \P\ — T on the graph starting from 
a vertex x — xo £ X and ending in v — xt € L such that if < j < T 
then \ L. Such a walk will be referred to as walk from x to L and 

denoted by P : x — > L or P : x — > v if we want to specify the endpoint. 
(Thus we allow the x € L case. It is implicitly understood in the notation 
that the walk does not visit L except for the last vertex and maybe the 
first.) For any such walk P we define a weight by the formula 

\p\-i 

S{P) = ^ d X] d Xj+1 (1 - A - QW(xj j) ' (14) 

Then we construct our extension of / using the following expression. If 
x G X \ L then put 

f( X )= J2 f(x\ P \)s(p) = ^2f(v) J2 s ( p )- ( 15 ) 

First of all, it is easy to see that this infinte sum will be absolutely 
convergent for large Q (and hence large A). Since W is strictly smaller 
than 1 outside of L, the terms 1/(1 — A — QW(xj)) will be uniformly 
bounded from above by some c/Q as Q goes to infinity. On the other 
hand, by induction on k we have 

fc-i 

(I l/\Jdx J d Xj+1 < Vd x d v . (16) 

P:x^v,\P\ = k 3 = 



G 



Hence if we sort the terms according to the length of P, we get a power 
series in 1/Q (or 1/A) with bounded coefficients. This easily implies ab- 
solute convergence. 

Let us consider now a walk P : x — > L and denote y = x\. Let R denote 
the walk y, x%, . . . , Xt, i- e. R is obtained from P by removing the first 
vertex, and P = xR. By definition s(R) = s(P)v/ 'd x d y (l — A — QW(x)), 
and obviously \R\ = \P\ — 1. Thus we can partition all paths P : x —> L 
according to their second vertex (denoted by x\) to arrive at the following 
identity: 

/(*) = E E /o=i*i+iM*fl) = 

xyGE R:y^L 

l 



E -7= E f(xm+MR) 



i-x-QW(x)^ E ^u- yR:y ^ 

= I V m (17) 

l-X-QW(x)) £^ E ^d^ 1 ' 

Multiplying both sides by 1 — A — QW(x) we see that (Hf)(x) = Xf(x) 
hence / is behaves like an eigenfunction of the Hamiltonian for every 
vertex in X \ L. We could actually use the same formula to "extend" / 
to L, only it is already defined there. Hence the condition that / is really 
an eigenfunction for H is equivalent to saying that the "extension" agrees 
with the starting value. Hence for every v £ L we get an equation 

/(«)= y, /w«( p ) = E/w E s ( p )- ( 18 ) 

P:v^L w£L Pw—*w 

Let us denote by Z VW (X,Q) = X/ P . S (P) the sum of the weights 
for all paths from v to w, as a function of A and Q. These functions 
form an \L\ x \L\ matrix Z — (Z vw ) and the conditions on the values 
/\l can be written simply as /\l = Z ■ Hence we have reduced the 
problem of analyzing eigenfunctions of an n x n matrix to understanding 
eigenfunctions of an \L\ x \L\ matrix whose entries are functions. This 
isn't a huge gain in general, but for the particular cases we are interested 
in, it does help a lot, as we shall see. 



2.3 Symmetric double wells — the proof of Theo- 
rem 11.81 

To illustrate the usefulness of the matrix Z(X, Q) we investigate the case 
when W has 2 global maxima: x and y, and there is an involution of G 
that preserves W. In this case there is a bijection between paths from x to 
x and paths from y to y. Also there is a bijection between paths from x to 
y and paths from y to x. These bijections obviously preserve the weight 
s(P). Hence Z xx = Z yy and Z xy = Z yx . This way the eigenvectors of the 
matrix Z are automatically (1,1) and (1,-1) with respective eigenvalues 
Z xx ± Z xy . This means that the eigenvectors for H are exactly as we have 
presumed in Section 11.11 Hence there is perfect asymptotic tunneling, 
exactly as described there. To determine the tunneling time we need 
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to understand that for fixed Q, how are the solutions of Z xx (Xi,Q) + 
Z xy (Xi, Q) = 1 andZ xx (\ 2 , Q) — Z xy (\2, Q) = 1 relate to each other. Let 
us introduce the auxiliary variable h = (1 — X)/Q. Now the weight of each 
path can be rewritten as 

If l-i 

1 1 

s(p) = w\ n ^—{h-wj^) (19) 

Let us write s' P = s' P (h) = s(P)(h — 1)Q' P L Then, since we know that 
h — > 1 as Q — > oo and that < W(v) < 1 when v is not a well, we get that 
s' P tends to some constant as h — > 1. Now, after multiplying by ft — 1, the 
equations for the eigenvalues of Z being equal to 1 can be written as 

P: X — V X P: X ^ry 

Denoting the solutions of each version by hi, ha, taking the difference of 
the two equations and dividing by hi — hi we get 

v s' P {hi) - s' P {h 2 ) Q~ d ^ v S'pjh!) + s' P {h 2 ) 

2^ (h 1 -h 2 )Q\ p \ hi -h 2 J—' Q\ p \~ d ' ( ' 

Since each s' P is clearly differentiable and the length of any x — > x path 
is at least 2 by definition, the first sum on the right hand side goes to 
as Q — > oo. In the second sum all the s'p's converge to some constants, 
and hence all terms where \P\ > d(x,y) go to 0, while those finitely many 
terms where \P\ — d(x,y) converge to constants. Hence Q d ^ x ' v ^/(hi — h 2 ) 
has to converge to a constant too. Thus we get hi — h% ~ l/Q d ( x ' y > and 
hence Ai — A2 ~ 1/Q d( - X,y ^~ 1 . Combining this with the observations in 
Section [O] we get that the tunneling time is of order Q d ^ x ' v ' ) ~ 1 . □ 



3 Simple potential 

In this section we are going to analyze the two and three-well scenarios 
where the potential wells are "singular", in the sense that W(v) = if v 
is not one of the wells and W(v) — 1 at each well. In this case the entries 
of matrix Z have particularly simple dependence on Q which makes the 
analysis possible. First of all, the weight function on paths becomes 

• ( ^^]STOr (22) 

when the starting vertex of P is not a well. For paths starting from the 
wells we get 

1-A 1 |PH1 1 

s(p) = i-A-Q ( i-A)i'i n Vdxjdxj+1 - (23) 
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This finally gives 

z vw (a, q) = y^t^q £ (t^tf £ n^4=- (24) 

|P| = fc 

The last part of the formula is independent of A and Q- Hence we can 
introduce the notation 

*«.(*)= e n /d ^ ( 25 ) 



LP|=* 



so we have 



Claim 3.1. For any a; 6 X t/ie uaZite of P xx (k) — PR(x,k) is exactly 
the probablility that the simple random walk started from x returns after 
k steps. 

Proof. For any closed path P : x = xo, xi, . . ■ , xt = x since the starting 
and endpoint are the same, we actually have 

T-l T-l 

n 7 ===nd L - (27) 

and this is exactly the probability of the random walk traversing P. Hence 
the sum for all P of length k is exactly the probability of return in k 
steps. □ 

We are interested in when 1 is an eigenvalue of the matrix Z(X, Q) = 
(Z VW (\,Q)), but this is clearly equivalent to (1 — A — Q)/(l — A) = 1 — 
Q/(l — A) being an eigenvalue of the transformed matix Z(X) which consist 
of entries 

(^ = £(T^P (28) 

Summarizing what we have computed until now: if Q is large, then the 
Hamiltonian will have two eigenvalues close to — Q. The corresponding 
eigenfunctions restricted to the wells will be eigenfunctions of the matrix 
Z with eigenvalue 1 — Q/(l — A). On the other hand if you start with 
an eigenfunction of Z with eigenvalue 1 — Q/(l — A) then extending it by 
formula (|15l) we get an eigenfunction of H with eigenvalue A. We are in 
very good shape now, as for any eigenfunction of Z(X) there is a unique 
Q that makes it into an eigenfunction with eigenvalue 1 — Q/(l — A) and 
the preceeding argument ensures that Q w —A if A is large enough. This 
way we can get rid of the Q parameter and only work with Z{X). This 
matrix is symmetric, since for any path P : v — » w the reverse P is a path 
from w to v with the same length and weight (the latter is so because we 
are using the symmetrized Laplacian). 
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3.1 Double wells - the proof of Theorem 11.31 

Let us start with the classical case, when there are two wells, x and y and 
the matrix Z(X) is 2 x 2. The first question we have to answer is whether 
there is tunneling or not. This depends on how the two eigenvectors of 
Z(X) behave. Let us denote the two coordinates of an eigenvector by 
f\{x) and f\(y), noting the dependence on A. Since Z xy (X) 7^ 0, neither 
coordinates of the eigenvector will be zero. Thus in both equations coming 
from (|18[> we can divide by term on the left hand side and combine the 
two equations into one, obtaining 

Z m (X) + fy^\z xy {\) = fy&Z yx (X) + Z yy {X). (29) 

a 0*0 fxw 

Using the symmetry of Z and putting a = a(f) — f\(x)/f\(y) we get 

Z X y ( OL ] = Z XX Zyy . (30) 

V a J 

Since our eigenvectors are normalized, we are only interested in the 
limiting behavior of a as A — > —00. Since there are two eigenvectors and 
they are orthogonal, if the other eigenvector is g then a(g) = —l/a(f). 
Hence the limiting behavior of a(f) and a(g) together as an unordered 
pair is clearly determined by the behavior of a — 1/ot. Hence the only 
thing we have to understand is 

lim fa-lU hm ^W-^W. (31) 
\ a J Z xy (\) 



Further simplifying notation by writing t = 1/(1 — A) we have to analyze 
the limit 

The denominator's leading term is clearly t d ^ x,v ' where d(x,y) denotes 
the distance between x and y in the graph. Depending on the graph 
three different limiting behaviors are possible. If for some k < d(x, y) the 
difference P xx (k) — P yy (k) is non-zero (i. e. the cospectrality of the wells 
is less than d — 1), then the limit will be ±00, hence the two eigenvectors 
tend to (0, 1) and (1,0). In this case the initial state of the particle is 
almost an eigenstate itself, so the particle will remain very close to its 
initial state, there is no asymptotic tunneling at all. 

If for all k < d(x,y) the difference P xx (k) — P yy (k) is zero (i. e. when 
the wells are d-cospectral), then the limit is 0, hence the two eigenvectors 
tend to (1/V2, l/\/2) and (1/^2,-1/^2). This is exactly the setting 
of Section 1 1.1 1 As we have seen in Section 12.11 all but two Cfc's in (0 
converge to 0, and we have just shown that the remaining two converge 
to the same value 1 / v2- Hence according to the argument in Section 11.11 
there is asymptotically perfect between x and y. 

Finally in the case when the first non-zero difference P xx (k) — P yy (k) 
is for k = d(x,y), that is, the cospectrality of the wells is exactly 
d — 1, the limit will be a non-zero constant c = (P xx (d(x,y)) — 
P yy (d(x,y)))/P xy (d(x,y)). Hence the two eigenvectors will converge to 
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(a, b) and (—6, a) where a/b — b/a = c,a 2 + b 2 — 1. In this case a ^ b and 
in ((7|) the two non-vanishing coefficients stablize to a and —6. Hence we 
get 

I \ ^ it\ 2 ( it(\i—\ 2 ) 2 | j,2\ 

ipt(x) ^ e I e k 1 -"a + b j . 

So maximal tunneling still occurs at time t — 7r/(Ai — A2) however it is not 
perfect tunneling. The probability of the state x never goes asymptotically 
below b 2 , and the probability of the state y never goes above a 2 — 1 — b 2 . 

To finish the proof of Theorem 11.31 we have to estimate the tunneling 
time. From the arguments in Section 11.11 it is clear that the tunneling 
time is approximately 7r / 1 Ai — A2I where Ai and A2 are the eigenvalues 
corresponding to the two eigenfunctions concentrated on the wells. The 
difficulty is that until this point we "fixed" A, chose an eigenvector for 
Z(X). We didn't have to worry about the eigenvalue, because we could 
always find a suitable Q for which everything worked in the end. Of course 
since Z has two eigenvalues, so we actually get two choices for Q each of 
which gives one of the eigenvectors in return. Now we have to reverse 
things, and compute the two possible A's as a function of Q. To do this, 
first we describe how the two possible Q's behave as the function of A. 

Let 1/1,2 denote the two eigenvalues of Z. Then, since v = 1 — Q/(l — A) 
has to hold, the possible choices for Q are given by Qi,2 = (1 — A)(1 — 2/1,2). 
The i/1,2 are given as the solution of a quadratic polyomial, so we can write 
them down explicitly: 

Zxx ~\~ Zyy zil \l (Z X x Zyy) 2 ~~\~ i*2Zxy) 2 

"1,2 = 1 2 , (33) 

hence both roots are real and we can compute: 

h{\) = Qi - Q 2 = (1 - A) fy/ '(Zxx - Zyy)' + (2Z xy ) 2 ^ . (34) 

We have already seen that Q\ is a power series of the form Q\ (A) = A+co + 
ci/A+C2/A 2 + . • • , and by the definition of h we have Q2(A) = Qi(A) — h(\). 
Let us denote the inverses of Qi(A) and Q(A) by Ai(Q) and \2(Q)- We 
will see that h(X) is a power series in 1/A with no constant term and then 
by the following lemma we get that A2(Q) — Ai(Q) = h(Q) + o(h(Q)), so 
the tunneling time turns out to be n/(h(Q) + o(h(Q))). 

Lemma 3.2. Let f(t) = t + c + ci/t + c 2 /t 2 + . . . and let g(t) = fit) -hit) 
where h(t) is a power series in 1/t with no constant term. Then g^ x (t) = 
f-\t) + h{t) + o{h{t)). 

Proof. g~\t) = r l {f{g-\t))) = r X <J)(a-\t)) + h(g-\t))) = r\t + 
Kg-Ht))) = f- 1 (t)+f- 1 '(t)h(g- 1 (t)) + o(h(g- 1 (t))). But f- 1 ' (t) = 
1 + o(l) and g- l (t) = t + 0(1) hence hig- 1 ^)) = h(t + 0(1)) = hit) + 
o(h(t)). Putting these together we see that indeed g l (t) — f~ (t) + 
h(t) + o{h{t)). □ 

It can be seen from (134[) that the order of magnitude of h(X) is the 
larger of the magnitude of Z xx — Z yy and of Z xy . The latter is always 
l/\ d{x > v \ while the former is by definition 1/A°°^'^ +1 . There is tunneling 
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only if co(x, y) + 1 > d hence the magnitude of h is always governed by 

the 1/A d term, hence the magnitude of the tunneling time is l/h(Q) « 
Qd(x,y)-i ag claimed n 

3.2 Triple wells 

Let now G be a vertex-transitive graph and let the three wells be denoted 
by x, y, z. This time we have to analyze the eigenfunctions of a 3-by-3 
matrix Z(\) instead of a 2-by-2. Let us denote the pairwise distances 
between the wells by a — d(x,y) < b = d(x,z) < c — d(y,z). Henceforth 
we shall assume that the wells are pairwise 2a-cospectral. Though the 
general case could be done along the same lines, we confine ourselves to 
the study of sufficiently cospectral wells mainly in order to reduce the 
number of cases to check. Already in this special case we are able to 
exhibit interesting phenomena. As we shall see, there are different type 
of behaviors depending on the distances between the wells: (i) a < b < c, 
(ii) a — b < c, (iii) a = b = c. Before looking at the particular cases, we 
can further simplfy the matrix Z a little. 

If we subtract Z XX (X) from each diagonal entry in Z then we only 
change the eigenvalues, but not the eigenvectors of the matrix. By the 
cospectrality assumption Z yy — Z xx = 0(1/(1 — A) a ) — Z zz — Z xx . The 
modified matrix for the triple-well case looks like this (remembering that 
a<b<c): 

-, / / (A) /U„(A) \ 

Z-Z XX I=- — ■ /o(A) / 0+0 i(A) /c-a(A) (35) 

U A) V A-a(A) /c-a(A) / 0+0 »(A) J 

where a', a" > and each entry denotes a power series, the index in- 
dicating the first term that may have non-zero coefficient: f m (t) = 
X]fcL m c fc/(l ~ A) fc , such that c m =fc 0. By abuse of notation, even if 
there are equalities among the indices, the corresponding power series are 
allowed to be different. For large A the eigenvalues and eigenvectors of 
this matrix can be approximated by taking the limit of the matrix on the 
right hand side. Let 

/ /o(A) /U„(A) \ 

M= lim /o(A) /„+„/(A) /c-a(A) (36) 

V A-a(A) fc-a(X) f a+ a»W J 

denote the limit of this matrix as A — > oo. Note that M depends only on 
the graph and the position of the wells. Let ^i,/U2,/X3 denote the three 
eigenvalues of M and ipi, ip2,ip3 the corresponding eigenvectors. 

Lemma 3.3. Assume that M has three distinct eigenvalues. Then a 
neccesary condition for perfect asymptotic tunneling from y to z is that 
ipz{x) = Vs(y) + ~4>3,{z) = = ipi(y) - ipi(z) = ip2(y) - ip2(z) for some 
permutation of the eigenvectors. 

Proof. Assume that there is perfect asymptotic tunneling from y to z. 
Since M has no multiple eigenvalues, the three relevant eigenvectors of 
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the Hamiltonian restricted to the wells converge to the three well-defined 
eigendirections of M. Using the definition of perfect tunneling and com- 
pactness it is easy to see that we get three real numbers r\,r-2,rz and 
three unit complex numbers pi,p2,p3 such that (0,1,0) = Y2 r i' i l > i an( A 
(0,0,1) >>,'•.<.',• 

If either of the rj 's would be 0, then (assuming n is zero) we get that 
(0, p2, —1) = (p2 — P3)r3ip3. Since ip3 is real, this can only be if p2 = ±1. 
Since obviously p3 / p2 we may assume wlog that p% =fc — 1, hence pi = 1 
and ipz = c(0, 1, —1)- This proves the part of the statment involving ip3, 
the rest follows by orthogonality of the eigenvectors. 

Now we assume neither of the r/s is 0, and denote = rjipj. Next 
we show that 9j(x) has to be for at least one j. Assume this is not 
true. Then ^2^i( x ) = an d "12 Pj^ j( x ) = 0. From this we get that 
(pi — p2)^2(x) + (pi — P3)^3(x) = 0. Since none of the ^j(a;)'s are zero, 
this implies that pi is a convex combination of p2 and p3. For unit complex 
numbers this can only be if they are all equal - a contradiction. 

So we may assume ^3(3;) = ip3(x) — 0. But then ^i(x) = — ^2(2^) 7^ 
and pi"i/i(x) — —p2^2(x), hence pi — p2- Looking at y we get that 
Vi(y) + *2(y) + *a(y) = 1 while pi(*i(y) + * 2 (y)) + P3*3(y) = 0. 
This can only happen if pi = — p3 = ±1. We may again wlog assume 
that pi = 1 and p3 — —1. Finally taking the difference (0, 1 — 1) = 
(0, 1, 0) - (0, 0, 1) = 2r 3 ^3 we get that ^3(2) + Vfe(v) = 0. The rest of the 
statement follows from orthogonality again. □ 

Whether the above condition of is also sufficient for perfect tunneling 
is described in the next lemma. 

Lemma 3.4. Assume M has three distinct eigenvalues and that tl>3(x) = 
ip3(y) + ip3(z) — and that neither if>i(x) nor ip2(x) is 0. Let 

7=^^. (37) 

pi — p.2 

a ) If J — P/l * s rational where p £ 7L is odd while q £ Z is even, then 
there is perfect tunneling from y to z in time 0(Q a ~ 1 ). 

b) If j is irrational, then there is perfect tunneling from y to z, but the 
tunneling time f has to satisfy Q a_1 — o(f(Q)). 

c) Otherwise we cannot say for sure that there is perfect tunneling from 
y to z. However if there is, the tunneling time f has to be satisfy 
Q a = 0(f(Q)). 

Proof. Let v be an eigenvalue of Z with eigenvector ip. Let us recall that 
ip is the restriction of an eigenfunction of the Hamiltonian to the wells 
with eigenvalue A if and only if Q = (1 — A)(l — v). If we denote by 
Vj = Vj(X) : j = 1,2, 3 the eigenvalues of the matrix on the right hand 
side of (|35f) then on one hand pj = lim,\->cx) Vj, on the other hand the 
eigenvalues of Z are Vj/{1 — A) a + Z XX (X) and the solutions for Q are 
given by 

Qj = (1- A)(l- Z xx (\)) - 13/(1 - A)"" 1 : j = 1,2,3. (38) 
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Denoting the inverse functions by Aj(Q) : j = 1,2, 3 we see that the three 
relevant eigenvalues of the Hamiltonian are precisely the Xj(Q)'s. 

Since ^i(x) and ij)2{x) are non-zero, all three eigenvectors participate 
in the decomposition of the pure state y. Then from the proof of the 
previous lemma we can see that perfect tunneling takes places from y to 
z in time f(Q) if and only if for large Q there is a t £ [0, f(Q)] and a unit 
complex number p such that e" Al(Q) « e ,; * A2<Q) « p and e" A3<Q) « -p. 
(Here Aa(Q) corresponds to the eigenvector ^3 in the previous proof.) 
More precisely when 

limsup inf \ e ^Q)-^(Q))_ ll 2 + le u { x l{ Q ) -x 3(Q)) + ll 2 =0 {w) 
Q^oo te[o,/(Q)]' 

Let us introduce the lattice C = {(p, 9) £ Z 2 : p = 0(2), g = 1(2)} and 
denote by 

r = {(f MQ)-W) >t A 1 (Q)-A,(Q) )eRa;t6[0>/(q)]} (40) 

7T 7T 

the relevant curve traced out in R as t runs over the interval [0, f(Q)]- 
Then (|39|) is equivalent to 

limsupdist(£,r) = 0. (41) 

Q — >oo 

Let us analyze what happens to the curve r as Q — > 00. Using (|38|) 
together with Lemma 13.21 we can write 

A* (0) - (0) = (1 + 0(1/Q))K (0) - ifc (Q))/(l - Q) a -\ (42) 

If /(Q) — °{Q a ) then clearly T converges uniformly to a segment through 
the origin whose slope is 7 and whose length is of magnitude ,f(Q)/Q a ~ ■ 
First of all, if f(Q) = o{Q a ' 1 ) then T converges uniformly to a single 
point, namely the origin. As (0,0) ^ C, in this case there can not be 
perfect tunneling. 

a) If 7 = q/p is rational where q £ Z is odd while p £ Z is even, then 
the line with slope 7 passes through a point of C. Hence there is a 
constant K such that if f(Q) = KQ a ~ l then the limit segment of 
r will pass through C hence (|41|l will hold. Thus there is perfect 
tunneling in time 0(Q a ~ 1 ). 

b) If 7 is irrational, then the line of slope 7 will pass arbitrarily close to 
L. If f(Q) = 0(Q a ~ 1 ) then the length of the limiting segment will be 
bounded hence the lim sup in (|4ip will be strictly positive and there 
will not be perfect tunneling. On the other hand if Q a_1 = o(f(Qj) 
then F will arbitrarily approximate any finite piece of the limiting 
line, hence there will be perfect asymptotic tunneling. 

c) Finally, if 7 is rational but not odd-over-even, then the lattice will 
be separated from the limiting line. Hence if f(Q) = o(Q a ) then the 
limsup in (|41|l will be positive again. Thus tunneling cannot take 
place within time o(Q a ). 

a 
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Remark 3.5. The only place where we use that M has no multiple 
eigenvalues is that it implies that the eigenfunctions converge to the three 
eigenvectors of M. If M does have multiple eigenvalues but we can com- 
pute the eigenvectors of Z, and denote their limits by ^1,2,3, then the 
statement of the Lemma remains true, as the same argument works. 

It is not hard to see, that if (Ai(Q) — A3(Q))/(Ai(<2) — A2(Q)) is not a 
constant as a function of Q, then even in the last case there will always be 
perfect tunneling and the tunneling time will depend on the degree of the 
first non-zero term in (Ai(Q) — \s(Q))/ (Ai(Q) — ^2(Q)) — 7- On the other 
hand if this is a constant function and 7 is rational but not odd-over-even, 
then there is no perfect tunneling at all. 

Proof of Theorem \1.4\ We have to consider two cases: 
a < b < c : In this case 

/ ci \ 
M = ci , (43) 
\ / 

hence the eigenvectors of Z converge to (1, 1, 0)/v2; (1, — 1, 0)/\/2; (0, 0, 1) 
with respective eigenvalues ci, — ci, 0. This of course means that there is 
complete asymptotic tunneling between the two nearest wells, and there is 
no tunneling to or from the third well. Further we get that the tunneling 
time is asymptotically 7r/|Ai(Q) — A 2 (Q) = 2ci7r/(l — Q) a_1 . 



a = b < c : In this case 

M=| ci I . (44) 
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Let us denote r = \J c\ + c| . Then M's eigenvectors are 



^ = (45) 
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(0,-c 2 ,ci) 



The corresponding eigenvalues are —r, r and 0. Some computation from 
(|35[) shows that ^i(A) ~ 1/(1 — A) c ~ a . So if the particle starts from state 
x then, since (1,0,0) is a multiple of (t[>2 — ipa), we get that the particle 
tunnels to the mixed state in which the probability of being in y and z is 
both 1/2. 

The situation is different if we start the particle from y. By Lemma[3]3] 
it follows that if Ci 7^ C2 then there cannot be complete tunneling from y 
to z. So let us assume ci =01. Since r / 0, all conditions of Lemma 13.41 
are satisfied and I7I = 1/2, hence part a) applies and we get that there is 
perfect tunneling in time 0(Q a_1 ). □ 
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3.3 Equilateral triple wells 



We are still using all the above notations, in particular we are analyzing 
Z using (|35|) . Now we have 

(0 Cxy C X z \ 
Cxy Cy Z , (46) 
Cxz Cyz / 

where the c's are defined according to ([5]) and depend on the particular 
geodesies connecting the wells. 

Proof of Theorem \1.5\ All of the c's are strictly positive, hence M has no 
eigenvectors with a single non-zero entry. This implies that some partial 
tunneling always happens. The question is how to determine when perfect 
tunneling happens. It is easy to see that M is non-singular, hence its 
eigenvalues are non-zero reals. The characteristic polynomial is x 3 — (c xy + 
c xz + c yz) x ~ 2c X yCxzCyz- If this has multiple roots then it has a common 
root with its derivative and just looking at the signs, this common root 
can only be x — — (c xy + c xz + q; z )/3. But that implies <yc xy c X zC yz = 
\J c xy + c| z + qJ z /3, which can only happen if c xy = c xz = c yz since the 
c's are non-negative. 

Let us first assume this is not the case. Then there are no multiple 
eigenvalues hence there are three distinct eigenvectors, which are the limits 
of the corresponding eigenvectors of the right hand side of (|35[) . As before, 
we denote the three eigenvectors by ip\, ip2, ip3- 

Clearly, the condition of Lemma 13.31 for 4>3 is equivalent to c xy = c xz . 
This proves part a) of the theorem. It is also not difficult to compute that 

c - ± /c 2 +8c 2 

the eigenvalues are = — — ~>M3 = ~ c yz- Thus in this case 

_ 3c y z + \J c y z + 8c| z 
c yz ~\~ 8c| z 

Since the c's are rational numbers by construction, 7 is rational if and 
only if -y/ c yz + 8cl z is rational. 

Claim 3.6. If j is rational it cannot be odd- over- even. 

Proof. By multiplying through with the appropriate integer, we may 
assume that Cj^c^ are coprime integers. If c yz is odd then so is 
\J c yz + 8c| z , hence the numerator is even and the denominator is twice 
an odd number, hence 7 cannot be odd-over-even. So we may assume 
c y z = 2p, hence c xz is odd (as they are coprime). But then 2-y/p 2 + 2ca; Z 
cannot be an integer, since p 2 + 2c xz is 2 or 3 modulo 4, which means it 
cannot be a full square. This is a contradiction hence the claim is true. □ 

Lemma [3. 41 now implies parts b) and c) of the theorem, except for the 
case when c xy = c xz = c yz . In this last case the difficulty is that the 
eigenvectors of M axe not well-defined, hence the arguments where we 
use that the eigenvectors of Z converge to the eigenvectors of M do not 
work automatically. This can be overcome by choosing a sequence of Q's 
along which the eigenvectors of Z do converge, and denoting their limits 
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(which will of course still be pairwise orthogonal eigenvectors of M) by 
ipi,ip2>ip3- Now by the arguments of Lemma 13.31 we still get that if>3 has 
to be parallel to (0, 1, —1). Also we know that tfji is parallel to (1, 1, 1) 
since this is a 1-dimensional eigenspace of M. And then ^2 is parallel 
to (—2,1,1). Now the corresponding eigenvalues are 2,-1,-1 and the 
arguments of Lemma [3.4| go through word by word. As 7 = 1 in our case, 
and 1 is not odd-over-even, we get that the tunneling time has to satisfy 
Q a = 0(f(Q)) as before. This completes the proof of part c). □ 

Remark 3.7. Now we see that by Remark l3.5l the only way Coniecture ll.6l 
could be false is if there existed a graph with three of its vertices x, y, z 
forming an equilateral triangle such that c xy — c xz 7^ c yz and w 'c yz + 8c xz 
is rational, and at the same time the ratio (Xi(Q)— \$(Q))/ (Ai(Q)— Aa(Q)) 
is independent of Q. Though currently we are unable to prove it, we 
believe such graphs do not exist. 

3.4 Instability 

Lemma 3.8. // there is an order three symmetry of the graph that per- 
mutes the wells (in particular c xy — c xz — c yz ) then there is no perfect 
tunneling from y to z, in fact TC(y, z) = 4/9. 

Proof. Using the order three symmetry of the graph it is easy to see that 
Z X x = Zyy — Z zz and Z xy — Z yz — Z zx . Hence (1, 1, 1) is an eigenvector 
of Z, and the whole subspace orthogonal to (1, 1, 1) is a 2-dimensional 
eigenspace. This way in the decomposition of (1,0,0) there will only be 
two terms: 1/3(1, 1, 1) + 1/3(2, —1, —1). As before, these two eigenvectors 
are rotating at different speeds, so after certain amount time (not hard 
to see that the order of magnitude is Q° _1 ) their phases will be opposite, 
and the state of the particle becomes l/3(— 1,2,2), which will be the 
extremally tunneled state. Hence TC(y,z) = 4/9. It is also obvious that 
at any moment during the whole evolution the states y and z will be 
equally likely. □ 

Proof of Theorem \1. 7| Let us consider the following simple construction. 
Let our graph consist of three paths joined at a common vertex O. (So O 
has degree 3, the endpoints of the paths have degree 1 and all the other 
vertices have degree 2.) Let the three wells be the three neighbors of O. 
Let the length of the paths going through x, y, z be denoted respectively 
by a,6,c. If a = b = c then by the previous lemma we get that for any 
two wells the tunneling coefficient is 4/9. 

On the other hand if a 7^ b — c then there is perfect tunneling from y 
to z. To see this, first we have to understand the eigenvectors of Z Let us 
assume a < b, though the a > b case would work just as well. It is easy 
to see that Z xy — Z xz — Z yz and by symmetry Z yy = Z zz . Since up to 
length 2a — 1 the closed paths from all three wells are identical, but the 
paths of length 2a are not, 7^ Z xx — Z yy l/\ 2a . Thus Z is of the form 




(48) 
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where p 7^ r. Easy computation shows that the eigenfunctions of this 
matrix are (ci, 1, 1); (02, 1, 1); (0, 1, —1) and ci / / C2. It is also not 
hard to compute that all three eigenvalues are distinct. The eigenvalues 
of the corresponding M matrix are 1,-1,-1 and hence 7 = 1. Thus, 
by Remark 13.51 in order to be able to apply Lemma 13.41 and conclude 
that there is indeed perfect tunneling from y to 2, we have to check that 
(Ai — A3)/(Ai — A2) is not a constant function of Q. If it were constant, 
it would be equal to its limit as Q — > 00, which is precisely 7. But since 
7 = 1, this meant that A2 = A3 for all Q, which means that the inverse 
functions are also equal, and hence for the eigenvalues of Z we get that 
vi (A) = v$ (A). But this contradicts our previous observations. Hence 
there is perfect tunneling from y to z. □ 
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